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Abstract 

In the PHENIX experiment at RHIC, the electro-magnetic calorimeter plays an 
important role in both the heavy-ion and spin physics programs for which it was 
designed. In order to measure its performance in the energy range up to 80GeV, a 
beam test was performed at the CERN-SPS H6 beam line. We describe the beam 
test and present results on calorimeter performance with pion and electron beams. 
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1 Introduction 



The PHENIX experiment at RHIC started data taking in 2000, with relativis- 
tic heavy-ion collisions. Subsequent runs will also use polarized proton beams 
to carry out a program of spin physics. In PHENIX, the electro-magnetic (EM) 
calorimeter is the primary tool for measuring photons and electrons/positrons. 
In order to cover topics in both physics programs, e.g. a thermal photon mea- 
surement in relativistic heavy-ion collisions, and prompt photon, 7r° and weak 
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boson measurements in spin physics, the EM calorimeter needs to cover a wide 
energy range extending from a few hundred MeV to 80GeV. A goal of the spin 
physics program is to measure differential cross sections of prompt photon and 
7r°'s to an accuracy of 10%. A 2% accuracy in the calorimeter energy scale is 
required to achieve this for p T > lOGeV/c of interest, because the cross sec- 
tions fall steeply as their energy increases. The EM calorimeter was originally 
designed for relativistic heavy-ion physics. There are two kinds of calorimeter 
in the PHENIX detector. One is a Shashlik [1-3] type lead-scintillator sam- 
pling calorimeter (PbSc) and another is a lead glass calorimeter (PbGl) . Table 
1 shows their basic parameters. A "super-module" is composed of 12 x 12 chan- 
nels for the PbSc calorimeter and 4x6 channels for the PbGl calorimeter. The 
total EM calorimeter system in the PHENIX detector consists of the PbSc 
super-modules and the PbGl super-modules. 





PbSc 


PbGl 


radiation length (Xq) [mm] 


21 


29 


Moliere radius [mm] 


~30 


37 


channel 






cross section [mm 2 ] 


52.5x52.5 


40x40 


depth [mm] 


375 


400 


[*o] 


18 


14 


7] coverage 


0.011 


0.008 


4> coverage 


0.011 


0.008 


super-module 






number of channels 


144 (12x12) 


24 (4x6) 


sector 






number of super-modules 


18 (3x6) 


192 (12x16) 


total system 






number of sectors 


6 


2 


number of channels 


15552 


9216 


rj coverage 


0.7 


0.7 


4> coverage 


90°+45° 


45° 



Table 1 

Basic parameters of two kinds of PHENIX EM calorimeter. 

The calorimeter's energy resolution, linearity and hadron rejection had already 
been measured at BNL-AGS in the energy range up to 7GeV [4]. In order to 
extend these measurements to the energy range up to 80GeV, a beam test was 
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performed at the CERN-SPS H6 beam line in 1998. 

In this article, we describe the beam test and present results of the PbSc data 
analysis. 



2 Setup 

Figure 1 shows a setup of the beam test. One PbSc super-module and four 
PbGl super-modules were located at the H6 beam line, and both were tested 
with electron beams in the momentum range of lOGeV/c to 80GeV/c and 
7r + beams of 40GeV/c. Both kinds of calorimeter were placed on a movable 
platform to change positions and angles of the incident beam on the calorime- 
ter. Delay-line Wire Chambers (DWC) [5] were located just in front of the 
calorimeter for measurements of the vertical and horizontal beam incident 
position. Incident-position dependence of the energy deposit were measured 
and corrected using the DWC. Two scintillators (SI and S2) were used as 
trigger counters, and two other scintillators (muon counters) were set behind 
iron blocks to identify muons in the beam. There was a Cerenkov Differential 
counter with Achromatic Ring focus (CEDAR) further upstream of the SI for 
electron identification. 

We used the lOGeV/c muon beams for channel-by-channel gain adjustment 
of the PbSc super-module channels in addition to the electron beams. For 
time-dependent gain drift correction of the PbSc calorimeter, we used a laser 
monitoring system [6]. 

The DWC has good position resolution (0.2mm) and high single-particle detec- 
tion capability (2xl0 5 particles/sec). It consists of one anode-wire plane and 
two cathode wire planes which surround the anode plane. The cathode planes 
have 2mm wire spacing. Their wires are connected with a delay-line through 
which signals are read out by a TDC module. Timing information corresponds 
linearly to position information. The active area is 100x100 mm 2 , and the area 
with linear response is 80x80 mm 2 

A beam trigger was composed of a SI signal ("SI"), a S2 signal ("S2") and a 
coincident signal of two muon counters ( "//" ) . An electron trigger was made 
by ST <E> S2 ® p, and a muon trigger was made by SI <E> S2 <g> /x. Triggers for 
pedestal measurement and for the laser monitoring system were used to take 
those data between beam spills. 

We used two different HV settings; 

• a "normal" HV setting (1.23-1.29kV) for energy measurements up to 80GeV 



3 



Om 



-- S1 



CERN SPS 
H6 beam line 



5m 



S2 



- DWC 



10m 




PHENIX 
PbSc+PbGI 



movable platform 



15m 



T. 



iron block 



muon counters 



Fig. 1. Setup of the high energy beam test at CERN. 
• a "low" HV setting (1.13-1.19kV) for energy measurements up to 160GeV. 

To readout PMT signals from the calorimeter, we used front-end electronics 
(FEE) from the CERN experiment WA98[7]. 



3 Analysis 

The deposited energy in each channel was calculated by multiplying the ADC 
count by a calibration factor, C(t) (GeV/count). The calibration factor has 
time dependence. We parameterize the time dependence by an initial gain 
factor, G (GeV/count), and a gain drift, D(t). 



The time dependent factor, D(t), is defined to be 1 at the time of the muon 
calibration run. 

The ADC count is derived from low and high gain ADCs which are both im- 
plemented to cover a wider dynamic range. We found a 1% difference between 
the ratio of high-to-low gain in electron trigger events and in laser trigger 
events. To determine the channel dependent high-to-low gain ratio, we use 



C(t) = G x D(t). 
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the measured ratio in electron beam trigger events for the 40 channels which 
have hits in the electron beam data. For the other channels, we measure the 
ratio with the laser trigger events. These channels have only a small contri- 
bution to the electron beam energy measurement hence a small effect on the 
systematic error. The average value of the high-to-low ratio is 7.8. The ADC 
value is derived from the low gain ADC when the low gain ADC count is larger 
than 90 counts. We also found there is no time dependence of the high-to-low 
ratio during whole run. The systematic error caused by the high-to-low ratio 
is less than 1% for those channels which have the electron beam data. There 
is a negligible contribution of the other channels to the systematic error. 

The gain factor is adjusted channel-by-channel by using muon trigger events. 
In order to identify muons, the following selections are applied. 

(1) The channel which has the largest energy deposit in all channels must 
have more than 80% of the energy sum of all channels. 

(2) The number of channels which have energy deposit more than 130MeV 
must be zero or one. 

When a muon beam penetrates one channel longitudinally, the most probable 
energy deposit is about 300MeV. Selection 1 requires that there are some hits 
which make a peak on that channel. Selection 2 requires it is a minimum- 
ionizing single particle and rejects background from other kinds of particles, 
and multi-hit of particles. By requiring both selections, muons which penetrate 
one channel are selected. 

After the muon selection, we have more than 100 muons in each of the 40 
channels. We adjust gains of these 40 channels so that the MIP peak position 
is at the same energy. The peak position is determined to a precision of 2- 
3%. In order to improve precision of the gain adjustment and to obtain the 
gain factor for the other channels, we use electron beam trigger events. The 
remaining errors of channel-by-channel adjustment of the gain factor is 3% in 
total. These errors are statistical ones. The systematic errors are smaller than 
these. 

An absolute value of the gain factor to provide correct energy scale is fixed at 
the electron beam energy of 20GeV. The result shows the average value of the 
gain factor is 110 (count/GeV). 

The time dependence of the gain was obtained using laser trigger events. The 
time variation of the laser amplitude was less than 3% reflecting the stability 
of the laser output. To monitor the fluctuation of the laser output, we use a 
truncated mean of the laser amplitude. 

The gain drift factor obtained with this method works reliably over periods 
of order a few hours. Between some sets of runs the gain drift is normalized 
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by using the beam energy of 20GeV in the period. The accuracy of the beam 
energy is 1% at 15GeV[8]. 

In this analysis, the total deposited energy is defined as a sum of energies in 
the 5x5 channels centered on the channel with the maximum energy deposit. 
The total energy is corrected by a position dependent factor. 

The upper-left figure of Fig. 2 shows the hit position dependence of the en- 
ergy sum in 5 x 5 channels for the 20GeV electron beam. In this figure, the 
coordinate (X, Y) shows a hit position in one channel obtained by the DWC. 
The hit position (0, 0) presents the center of the channel and (1, 1) presents 
the edge of the channel. The position dependence is fitted as shown in the 
upper-right figure of Fig. 2 by the following formula. 

1 + a x (X 2 + Y 2 ) + b x (X 4 + Y A ) + c x X 2 ■ Y 2 (1) 

We obtained a best fit with the following parameters; 

a = -0.3079 
6 = 0.3643 
c= -0.02894 

We use these parameters to correct for the position dependence of the energy 
measurement. The lower-left figure shows the deviation of the energy sum from 
the fitted hit position dependence and the lower-right figure shows a projection 
of the deviation. The deviation is 0.5% of the energy sum. The systematic error 
remaining after the position dependence correction is evaluated to be 0.5%. 



4 Results 

Figure 3 shows the efficiency for the 40GcV positron beam when we require a 
measured energy deposit greater than E cut , and the pion rejection power for 
the 40GeV ir + beam obtained with the same cut. 

The energy resolution at each energy point was obtained by fitting a Gaussian 
distribution within ±2er around the electron peak. The electron trigger events 
contained a 10% pion contamination. The contamination in the electron peak 
region is less than 0.1% because only 1% of the pions deposit more than 90% 
energy in the calorimeter. The x 2 of the fitting is reasonably small. 

Figure 4 shows energy resolution obtained by both beam tests at CERN and 
BNL. They can be fit with linear or quadratic expressions. Only statistical 
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Fig. 2. Upper-left: Hit position dependence of the 5x5 energy sum for the 20GeV 
electron beam, where (X,Y) shows a hit position obtained by the DWC. Up- 
per-right: Fitted hit position dependence. Lower-left: Deviation of the energy sum 
from the fitted position dependence Lower-right: Projection of the deviation. 

errors are taken into account in the fits. We estimate an additional 1% sys- 
tematic error based on the reproducibility of the measurements at each energy 
point. The results of the fits are 



a E /E = 1.2% + 
a E /E = 2.1%© 



6.2% 



y/E(GeV) 

8.1% 
y/E(GeV) 



where © denotes a square of the quadratic sum, a © (3 = ^/a 2 + f3 2 . They are 
valid in the energy region of 0.5GeV to 80GeV with 1% systematic uncertainty. 

Figure 5 shows the residual (measured energy with the calorimeter less the 
beam energy, divided by the beam energy) of the energy sum in 5 x 5 channels 
versus the beam energy. We see that the calorimeter response is linear within 
2% systematic uncertainties in the energy region from 20GeV to 80GeV. There 
is some indication of a 2% deviation from linearity at lOGeV, however this is 
within our systematic errors. Such a deviation cannot be due to one of the 
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Fig. 3. The efficiency for the 40GeV positron beam when we require a measured 
energy deposit greater than E cut , and the pion rejection power for the 40GeV ir + 
beam obtained with the same cut. 

following corrections; the gain drift correction, pedestal subtraction, high-to- 
low ratio correction, or the energy sum of 5 x 5 channels. We considered the 
following possible sources; 

• run-time problems in the monitoring system 

• linearity of the WA98 FEE 

• linearity of the PMT 

• an inherent non-linearity in the calorimeter due for example to the interplay 
between light attenuation in the wave-length-shifter (WLS) readout fibers 
and longitudinal shower leakage beyond the calorimeter. 

Linearity of the WA98 electronics which was used to digitize signals from the 
calorimeter was investigated. We found that it has linear response within 1%. 
Linearity of the PMT had also been investigated and confirmed to have linear 
response within 2% [9]. The light leakage in the WLS fiber is evaluated by 
the simulation program. The energy leakage is evaluated to be 1% at lOGeV 
and 4% at 80GeV [10,11]. These effects tend to cancel one another, so we do 
not expect the intrinsic non-linearity to be as large as 2%. In summary, the 
only place where the measured residual approaches the limit of 2% is at 10 
GeV. Nevertheless, the residual is consistent within our systematic error of 
being linear at this energy. We've examined a number of possible causes for a 



8 



14 
12 
10 
8 
6 
4 
2 



a BNL0.5-5.0GeV/ce 



CERN10-80GeV/ce 



\ 1 



1.2%+6.2%/%(GeV) 
2.1%e8.1%/NE(GeV) 



_l I I I I I 



_l I I l_ 



_l I I I I I I 



10 



GeV 



10' 



Fig. 4. Energy resolution obtained by both beam tests at BNL and CERN. A dashed 
line shows the result of fitting by a linear formula, ge/E = 1.2% + 6. 2%/ \J E(GeV). 
A dashed dotted line shows the result of fitting by a quadratic formula, 
o E /E = 2.1% © 8.1%/y/EjGeV). 

deviation from linearity but none has been found. 

Position resolution of the beam hit position is evaluated with the logarithmic 
method [12]. In the logarithmic method, the position is determined by the 
following formula. 



7~1 Cj x x i 



X 



i=\...N 



i=l...N 



where Xj denotes a center of each channel in the horizontal direction. Similarly, 
Y is denned in the vertical direction. The weights Cj are; 



R^Maxl^Ei/EtaoA 
Ci = Max[0,log(Ri)+Co] 

where E tota i is the total energy, E tota i = E^ and C is a constant. A larger 

i=l...N 

value of a and (3 is expressed by Max[a,0\. 
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Fig. 5. Linearity of the 5x5 energy sum for both beam tests at BNL (left) and CERN 
(right). Solid lines show total systematic uncertainties in the analysis. Absolute level 
of two beam test is not normalized. 

The deviation in a short period shows that systematic uncertainty in the 
logarithmic method is 2mm. 

Figure 6 shows the position resolution obtained by both beam tests at CERN 
and BNL. The points can be fitted by a formula; 

. 5.9(mm) 
a x (mm) = 1.4 mm H :. 
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We measured the energy resolution, linearity and position resolution of the 
PHENIX EM calorimeter in a test beam at CERN. For the PbSc calorimeter, 
we obtained energy resolution of; 

a E /E = 1.2% + 6 ' 2% 
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Fig. 6. Position resolution obtained by both beam tests at BNL and CERN. A 
dashed line shows the result of fitting, 1.4 mm + 5.9 mm / W E(GeV). 



a E /E = 2.1%© 



1.1% 



E{GeV) 



and position resolution; 

5.9(mm) 

<y x {mm) = 1.4(mmJ + 



E(GeV) 



A major purpose of the test was to investigate the performance of the calorime- 
ter in the energy range up to 80GeV and, in particular, the linearity of response 
versus beam energy. Since our goal in PHENIX is to measure prompt photon 
and 7T° production cross sections with the calorimeter within 10% errors it is 
important to understand the linearity of the calorimeter at the level of 2%. 
In the analysis of the PbSc calorimeter, the response is found to be linear to 
approximately this level. 
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